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ABSTRACT 

I investigate the origin of the observed correlation between a gamma-ray burst's vF^, spectral peak 
£'pk and its isotropic equivalent energy E{sa through the use of a population synthesis code to model the 
prompt gamma-ray emission from GRBs. By using prescriptions for the distribution of prompt spectral 
parameters as well as the population's luminosity function and co-moving rate density, I generate a 
simulated population of GRBs and examine how bursts of varying spectral properties and redshift 
would appear to a gamma-ray detector here on Earth. I find that a strong observed correlation can be 
produced between the source frame i?pk and Eiso for the detected population despite the existence of 
only a weak and broad correlation in the original simulated population. The energy dependance of a 
gamma-ray detector's flux-limited detection threshold acts to produce a correlation between the source 
frame iJpk and -Eiso for low luminosity GRBs, producing the left boundary of the observed correlation. 
Conversely, very luminous GRBs are found at higher redshifts than their low luminosity counterparts 
due to the standard Malquest bias, causing bursts in the low -Bpk, high E^^a regime to go undetected 
because their iJpk values would be redshifted to energies at which most gamma-ray detectors become 
less sensitive. I argue that it is this previously unexamined effect which produces the right boundary 
of the observed correlation. Therefore, the origin of the observed correlation is a complex combination 
of the instrument's detection threshold, the intrinsic cutoff in the GRB luminosity function, and the 
broad range of redshifts over which GRBs are detected. Although the GRB model presented here is a 
very simplified representation of the complex nature of GRBs, these simulations serve to demonstrate 
how selection effects caused by a combination of instrumental sensitivity and the cosmological nature 
of an astrophysical population can act to produce an artificially strong correlation between observed 
properties. 

Subject headings: gamma rays: bursts — galaxies: star formation 



1. INTRODUCTION 

The search for empirical correlations between observ- 
able parameters has long been an important path to 
understanding the underlying nature of astrophysical 
sources. Such attempts at correlating data, though, are 
inherently risky, as without careful consideration of ob- 
servational biases, one can be lead to false conclusions 
regarding the strength and nature of observed correla- 
tions. In this paper I consider the origin of the much 
discussed correlation between a gamma-ray burst's vF^ 
spectral peak (-Epk) and its isotropic e g uiva. lent energy 
{Eiso), first reported by |Amati et ah] ( |2002[ ). The im- 
portance of such a correlation cannot be understated, as 
a tight relationship between a GRB's spectral proper- 
ties and total energetics would allow for the distance of 
cosmological GRBs to be determined from gamma-ray 
data alone, opening the possibility of using high redshift 
GRBs for cosmological distance-scale applications. 

Discussions regarding the origin of the i?pk — i^iso cor- 
relation ha ve had a long and cont entio us history in the 



literature. Nakar & Piran (2005) and Band & Preece 



([2005) showed that a sigmticant traction of BAl'SE de 
tected GRBs without known redshift could not be con- 
sistent with the relation at any redshift, indicating that 
the true correlation may in fact be much broad er than 
the one originally found by Amati et al. (2002). Like- 
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wise, [Butler et ah (|2007[ |2009[ ) showed that hard and 
under- luminous GRBs detecte d by the Swift s pacecraft's 
Burst Alert Telescope (BAT) |Gehrels et a"L][2b04j ) were 
inconsistent with the pie-Swijt Epk — i/'iso correlation. 
Furthermore, the normalization of the Swift based corre- 
lation shifts towards the detection threshold of the BAT, 
suggesting that the strength of the correlation is largely 
governed by the sensitivity of the collecting instrument. 

Questions regarding the influence of detector thresh- 
olds on the i?pk — Eiso correlatio n were examined in de- 
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tail fo r a variety of ins tru ments by Ghirlanda et al. ( 2008 1 
and Nava et al 
claim that evii 
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detections thresh- 
olds are clearly present in their spectroscopic samples, 
but that the strength and shape of the resulting correla- 
tion is not govern ed by these effects. 

More recently, Butler et al. (2010) applied a multi- 
variant analysis to over 200 Swift detected GRBs in or- 
der to investigate their intrinsic energetics distribution 
and spectral parameters. The authors found evidence 
for a very broad correlation between the source frame 
-Epk and Ejso, whic h only takes the form observed by 



Amati et al. 



events, mate 



( 2002 1 through the non-detec tion of weak 
mng the conclusions drawn by [Nakar fc Pi-| 
,ran, ( ^2005| ; [Band fc Preece[ ([2005t . 

m this paper, 1 examine the nature of the Epi^ — E[so 
correlation through a new approach that utilizes a pop- 
ulation synthesis code to model the prompt gamma-ray 
emission from GRBs in order to examine how bursts of 
varying spectral properties and redshift would appear 
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to a gamma-ray detector here on E arth. This work is 



similar to the analysis performed by Nava et al] ( 2008 1 
in that I try to quantify the degree to which detector 
thresholds produce the observed — i?iso correlation, 
although I take the additional step of using prescriptions 
for the population's luminosity function and co-moving 
rate density to examine how the cosmological distribu- 
tion of these events biases the observed properties of the 
detected population. 

I find that a strong observed correlation can be pro- 
duced between the source frame -Epk and i?iso for the 
detected population despite only a week and broad cor- 
relation being present in the original simulated popula- 
tion. The origin of this observed correlation is a complex 
combination of the instrument's detection threshold, the 
intrinsic cutoff in the GRB luminosity function, and the 
broad range of redshifts over which GRBs are detected. 

Any energy dependance on the flux limited detection 
threshold of a gamma-ray detector acts to produce a cor- 
relation between the inferred source frame Epy^ and Eiso 
for low luminosity GRBs, producing the left boundary of 
the observed E'p^ — i?iso correlation. Conversely, very lu- 
minous GRBs are relatively rare and therefore probabil- 
ity dictates that we would have to look to higher redshifts 
in order to detect these events. Therefore, GRBs in the 
low E'pk, high Eiso regime would go undetected because 
their Epk values would be redshifted to energies at which 
most gamma-ray instruments become less sensitive, pro- 
ducing the right boundary of the observed correlation. 
The combination of these selection effects act to produce 
a relatively tight correlation between i?pk and Eiso hi the 
population detected by the instrument, despite the ex- 
istence of a much broader intrinsic correlation that goes 
unseen by the observer. 

I present an overview of the population synthesis code 
used in this paper in Section 2.1, followed by a more in- 
depth description of the code in Sections 2.2, 2.3, and 
2.4. I present the burst demographics of a simulation of 
40,000 GRBs in Section 3 and discuss the implications 
of these results in regards to the observed Ep^^ — Eiso 
correlation in Section 4. Throughout the paper, I will 
refer to the source and observer frame i^Fi, peak as -Bpk 
and i?pk,obs respectively. 

2. GRB MODEL 

2.1. Model Overview 

Gamma-ray burst continuum spectra can evolve quite 
dramatically over the course of a burst. This evolution 
is generally characterized by an overall softening of the 
spectra to lower energies, with the peak of the vF^ spec- 
trum evolving through the detector bandpass over the 
duration of the burst. This evolution will be delayed by 
the effects of time dilation for GRBs at high redshifts, re- 
sulting in a longer observed spectral lag between the high 
and low energy channels and a broadening of the pulse 
profile. At the same time, the observed GRB flux falls 
as a function of increasing luminosity distance. The net 
effect is that soft and faint emission will become increas- 
ingly difficult to observe with gamma-ray detectors such 
as Swift or Fermi 's Gamma-ray Burst Monitor (GBM) 
(Meegan et al. 2009), which suffer large drops in sensi- 
tivity at energies below ~20 KeV. The observed bursts 
properties are therefore a complex convolution of the ef- 



fects of cosmological redshift and detector sensitivity and 
hence I turn to simulations to obtain a better idea of how 
these bursts would appear in the observer frame. 

Two empirical correlations form the basis for the GRB 
model that I have developed to investigate this question. 
The first is the hardness-intensity correlation or HIC, 
which relates the instantaneous hardness of the spectra 
and the instantaneous energy flux Fe, within individual 
pulses. For the decay phase of a pulse, the most common 
behavior of the HIC is a power-law relationship between 
Fe and the peak of the vF^, spectrum, -Bpk.obs, of the 
form Fe cc E^^y_ ^j^^ , where r/ is the HIC power-law index. 
The second correlation is the hardness- fluence correlation 
or HFC (Liang & Kargatis 1996) which describes the 
observation that the instantaneous hardness, or i?pk,obs, 
of the spectra decays exponentially as a function of the 
time-integrated flux, or fluence, of the burst. The HFC 
can be stated as i?pk = i?oe~*/*", where is the 

photon fluence integrated from the start of the burst and 
$ is the exponen t ial dec ay constant. 



Kocevski et al. (20031 show that both the HIC and 
Hl^'C correlations can be produced through simple rela- 
tivistic kinematics when applied to a spherical shell ex- 
panding at relativistic velocity. The curvature of a rela- 
tivistic fireball would make the photons emitted off the 
line of sight delayed and affected by a varying Doppler 
boost due to the increasing angle at which the photons 
were emitted with respect to the observer. Simulating 
a GRB pulse using only these two correlations to de- 
scribe the e volution of a Band photon spectrum (Band 



et al.] 1993 1 as a function of time reproduces the I'ast 
Rise Exponential Decay (FRED) pulse shape that is so 
ubiquitous in GRB data. 

With this time-resolved spectral model, I can simulate 
a GRB at a variety of redshifts and quantify the effects 
that redshift have on the inferred pulse properties. In or- 
der to convert our modeled photon spectrum into a count 
spectrum and eventually a count light curve, I take the 
simulated photon spectrum and fold it through an instru- 
ment response matrix. I then add a Poisson distributed, 
energy dependent, background spectrum derived from 
the median backgrounds of a sample of detected bursts. 
From this count -|-background spectrum, I can then cal- 
culate the resulting count light curve that the instrument 
would have produced for a given input pulse model. 

Using the simulated count light curve and spectra, I 
can then back out the burst's trigger significance, dura- 
tion, hardness, peak flux, and time integrated spectral 
parameters as would be inferred by the observer. Using 
the known redshift, I can then use these properties to 
determine the inferred source frame Spk and obtain an 
estimate of Liso^ and i?iso- By generating a large number 
of GRBs following a realistic redshift distribution, I can 
examine the distribution of properties for the bursts that 
are detected by the instrument and compare them to the 
true source frame distributions known otherwise only to 
nature. 

Sections 2.2, 2.3, and 2.4 describe the model outlined 
above in more detail. The model is currently written in 
a combination of IDL and Python and is available upon 
request]^ 

^ http:/ /www. kocevski.com/GRBModel/ 
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Fig. 1. — The hard to soft spectral evolution of i?pij as a function 
of time as prescribed by relativistic kinematics for a spherical shell 
traveling towards the observer with high Lorentz factor. 

2.2. Modeling Individual FRED Pulses 

In order to model a single FRED pulse as it would 
appear in the source frame of the GRB, I start by as- 
suming an empirical Band photon spectrum and evolving 
its £^pk and peak energy luminosity Le at Ep]^ through 
time. Drawing upo n the extensive BAT SE spectral cat- 
alog presented by (Preece et al. 2000), I set the low 



and high energy power-law indices of the simulated spec- 
trum to match the median BATSE determined values 
of apk — —1.1 and /3pk = —2.3 respectively. Likewise, 
I draw our initial ii'pk.o from a log-normal distribution 
centered at £^pk,obs ~ 200 x (1 -I- z) keV and falling off 
sharply above 300 x (1 -I- z) keV. Here, z represents the 
redshift at which we expect to see the largest number 
of events, assumed here to be the redshift at which the 
cosmic star formation peaks. 

I then evolve the spectrum's ij^pk and Le in time fol- 
lowing the relations derived by Kocevski et al. (2003') for 
the effects of a spherical shell expanding at relativistic 
velocity, namely: 



Epk{t) = Ep^fiV — 



LE{t) = Le,oV^ 



E, 



pk,0 



(1+i/T-ang) 



Le,o 



(l+V^a„gF 



(1) 



(2) 



where i?pk,o and Lea represent the initial E'pk and Le 
values in the emitting surface's co-moving frame. T> is 
the angle dependent Lorentz boosting factor for transfor- 
mations from this co- moving frame to the GRB's source 
frame (i.e the rest frame of the progenitor), given by: 



r(i-/3Ai) 



(3) 



The combined effect of the angle dependent Lorentz 
boost factor and the additional time of flight of pho- 
tons emitted off the line of sight is the delayed arrival 
of soft photons. Even in the idealized case of a curved, 
relativistic emitting surface, which emits for an instan- 
taneous moment in time and at a single frequency -Epk.o 
with luminosity Le.o, the observed profile will be a broad 
pulse with harder emission arriving first followed by the 
arrival of softer and weaker emission. Equations [l] and 
[2] only describe the decay phase of the pulse, so I take 
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Fig. 2. — The resulting fast rise exponential decay (FRED) pulse 
shape generated by the model by integrating the bolometric GRB 
spectrum at each time bin. 

the additional step of introducing a power-law rise term 
to the analytic function describing peak energy flux such 
that LE{t) ~ LE{t) X {t/toY. Varying the to for this 
power law component can adjust the resulting time to 
peak flux accordingly, allowing for the production of a 
FRED pulses with a variety of durations. 

An example of the hard to soft evolution that results 
from Equations [T] and [2] can be seen in Figure [1] Both 
£'pk and Le decrease with time, with the red dotted line 
representing the time integrated spectrum produced by 
summing the individual time-resolved spectra. The green 
dotted lines represent the nominal energy window of the 
BATSE instrument. A 3-dinicnsional array is used to 
store the time-resolved photon spectra calculated at one 
second intervals, producing a photon data cube as a func- 
tion of energy, time, and photon flux. An example of a 
bolometric photon light curve, produced by integrating 
each time-resolved photon spectra over energy, can be 
see in Figure [2] The resulting light curve reproduces the 
ubiquitous fast rise exponential decay profiles observed 
in GRB time histories. 

2.3. Population Synthesis 

The cosmological distances at which GRBs occur play 
a non-trivial role in the resulting properties seen in the 
observer frame. As the redshift of a GRB increases, 
-Epk.obs is pushed closer to the lower bound of the de- 
tector's energy window. As this happens, more of the 
GRB's low energy spectrum goes undetected by the in- 
strument, reducing the observed flux and fluence beyond 
what is expected from simply increasing the GRB's lumi- 
nosity distance. Therefore, it is important to simulate a 
population of GRBs with a realistic redshift distribution 
in order to properly account for this effect. 

To do so, I turn to the luminosity function (j){z) and co- 
movin g rate density p{z) estimates presented in Butlerj 
et al. (2010), where the authors utilize a multi- variant 
analysis on over 200 Swift detected GRBs with known 
redshift to infer (f>{z) and p(z) after taking into consider- 
ation the Swift-BAT detection threshold. They find that 
0(z) decreases sharply above Lcut ~ 53, and find no sig- 
nificant evidence for luminosity evolution as a function 
of redshift. They also conclude that the GRB co-moving 
rate density follows the global star formation rate out to 
z ^ 2 — 3 and then flattens at higher redshift. Based on 
there results, I adopt the following analytic descriptions 
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Eiso = 4.1 e53 ergs 
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Fig. 3. — The resulting count light curve obtained by folding the 
time-resolved photon spectra through a detector response and inte- 
grating the resulting individual time-resolved count spectra. The 
burst fluence and energetics are based on the duration obtained 
from the count light curve 

for (pL and p{z) 



dN 
d log L 



(4) 



dN 
dz 



cx + < 0.97 cx [l + zY^; 



0.97 < z < zi cx (1 + > zi (5) 



where aL = —0.22 and 6_l = —2.89 and (g p, ffi, .i?2) 
(3.4, -0.3, -8) for zi = 4.5 ([Butler et al!][20l0|. 



2.4. Converting Photons to Counts 

Gamma-ray detectors trigger on counts accumulated 
by the instrument over a specified time interval. The 
triggering algorithms used to measure the significance of 
these accumulated counts can be quite complicate d, as 
is the case with the Swift-BAT ( [Gehrels et al.|2004[ ), but 
was rather straight fo rward for the BAi'SE instrument 



(Fishman et al. 1994]) onboard the Compton Gamma- 
"B 



ray Observatory, for BATSE, the instrument could trig- 
ger whenever there was an 5cr excess of counts accumu- 
lated over three possible time timescales (64 ms, 256 ms, 
or 1024ms) compared to a running time averaged back- 
ground computed every 17 seconds. 

Whether a GRB triggers the instrument depends not 
only on the intrinsic luminosity of the burst and its dis- 
tance from Earth, but also on the location of i?pk,obs in 
the detector's energy window. In addition, a burst's spec- 
tral and temporal properties are inferred by observers us- 
ing count data produced by the instrument, which like- 
wise depends on -Bpk,obs- Therefore, I convert the result- 
ing photon model into a count model in order to asses the 
detectability of the burst as well as to measure the re- 
sulting burst properties as they would have been inferred 
by the observer. 

In order to convert the observer frame photon model 
into a counts model, I take the simulated photon data 
cube and fold it through an instrument data response 
matrix (DRM). For the purposes of this analysis, I use a 
BATSE response file that was generated for a real burst 
which occurred nearly at zenith for one of BATSE's Large 
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Fig. 4. — i?iso data plotted versus redshift for our entire simu- 
lated sample. GRBs with a trigger significance of more than 5.5 
cr above background in the 50-300 KeV energy range are shown as 
crosses whereas events that would have went undetected are shown 
as points. A standard Malquest type bias produces the observed 
correlation between Eiso and redshift. 

Area Detectors (LAD) . The DRM describes the distribu- 
tion of counts over the instrument's energy channels due 
to the arrival of photon of a given energy. The result is 
a counts data cube as a function of channel energy, time, 
and count rate. 

I then add a Poisson distributed, energy dependent, 
background spectrum derived from the median back- 
grounds from a sample of BATSE detected bursts to each 
count spectrum. A count light curve with a realistic 
background spectrum can then be produced as a func- 
tion of time by integrating each individual time-resolved 
count spectra. The resulting burst duration is then cal- 
culated using a Bayesian block algorithm. An example 
of the resulting count light curve, for a GRB placed at 
z = 2, and the subsequent Bayesian block reconstruction 
(solid green lines) can be found in Figure [2] 

Although only a small fraction of the GRBs that form 
the basis for the E'pk — fi'iso correlation were detected 
by BATSE, the broad energy range (20-1800 KeV) cov- 
ered by the instrument and its relatively straight forward 
trigger criteria make it extremely easy to model for this 
analysis. Furthermore, for purposes of the analysis pre- 
sented here, I do not aim to reproduce the exact form of 
the observed Spk — i?iso correlation, but instead try to 
investigate how the combination of instrumental sensi- 
tivity and the cosmological nature of GRBs can work to 
produce an artificially strong correlation that does not 
exist in the source frame of the population. 

3. SIMULATED DATA SET AND ANALYSIS 

Using the GRB pulse model and population synthesis 
code described in Sections 2.2 through 2.4, I simulate 
a set of 40,000 GRBs placed at a variety of rcdshifts 
following the co-moving rate density p{z) described in 
Equation [5] For each burst, the simulation code gener- 
ates a time integrated PHA-I file by summing the counts 
data cube over a source duration as determined through 
the Bayesian block algorithm. These PHA-I files, along 
with the original DRM used by the model, can then be 
read into XSPEC to obtain time-integrated spectral fits 
and photon and energy flux estimates in the 25-300 keV 
energy range, F and Fe respectively. 

Using the XSPEC determined spectral parameters 
along with the measured energy flux and the Bayesian 
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Fig. 5. — The redshift distributions for all simulated (solid) and 
detected (dashed) bursts shown in comparison to the redshift dis- 
tribution of Swift detected bursts (dotted). 
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Fig. 6. — The distribution of Ejso for all simulated (solid line) and 
detected (dashed) bursts in comparison with the Swift determined 
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block determined duration, I can calculate the energy 
fluence Shoi integrated over an energy range of 10 to 
10000 keV. Using the recorded redshift, I can estimate 
the inferred isotropic equivalent energy £"130 and lumi- 
nosity Liso, k-corrected to a standard energy range of 10 
to 10000 keV, in the source frame. These inferred pa- 
rameters are in addition to the true bolometric Eiso and 
Liso, known to the simulation, but otherwise unknown to 
the observer. 

I find that although systematic differences do exist be- 
tween these inferred energetics estimates and their true 
values due to instrumental effects, these differences do 
not heavily influence the results presented below. I re- 
serve a more detailed discussion of the systematic biases 
between a burst's inferred and true properties that are 
introduced by such instrumental biases for a future pa- 
per. In this work, I will focus solely on the duration, 
5'boi, £'pk, Eiso and Ljso inferred by the user from the 
counts data for bursts that were detected by the instru- 
ment and use the true values to display the properties 
of bursts that were not detected by the instrument and 
hence for which no inferred values could be measured. 

4. RESULTS 

The resulting Eiso data plotted versus redshift for our 
entire simulated sample can be seen in Figure |4j GRBs 
with a count rate yielding a trigger significance of more 



All Simulated GRBs 
Detected Simulated GRBs 




Epk (keV) 

Fig. 7.— The 

^pk,obs distribution for all simulated bursts (solid) 
and detected (dashed) bursts. The low end of the detected -Epk.obs 
distribution is influenced by the detector's limited sensitivity below 
~ 20 keV, whereas the high end of the detected -Epk,obs distribution 
reflects the true cut-off in the simulated population. 

than 5.5 a above background in the 50-300 KeV energy 
range are shown as crosses whereas events that would 
have went undetected are shown as points. A standard 
Malquest type bias produces the observed correlation be- 
tween Eiso and redshift, emphasizing the fact that lumi- 
nous bursts are statistically found at higher redshifts. 
The redshift distributions for all simulated bursts (solid) 
and detected bursts (dashed) are shown in figures pj with 
the redshift distribution of Swift detected bursts (dotted) 
also shown for comparison. The decline in the observed 
GRB population at high redshift, compared to the simu- 
lated sample, is largely an effect of the sensitivity of the 
detector. I find that without invoking luminosity evo- 
lution, the relative number of GRBs occurring at high 
redshift must increase in order to explain the number of 
high redshift detections. This qualitative observ ation is 
in agr eement with similar conclusions reported by [Butler] 



et al. 



( 12010) . 



I'he distribution of Eiso for all simulated events (solid 
line) is shown in Figure |6] along with the distribution of 
Eiso for the bursts that would have triggered the detec- 
tor (dashed line) . The simulated population ranges from 
10'''' < Eiso < 10^^, but the resulting observed popula- 
tion drops off below Eiso < 10^^ erg due to the detector 
threshold. Conversely, the high end of the detected Eiso 
distribution is not influenced by the instrument, and in- 
stead reflects the cut-off that I assume for the GRB lu- 
minosity function. The dotted line shows the Swift Eiso 
distribution as determined by Butler et al. (2010) for 
comparison. Although slightly shitted to higher values, 
the Swift observed i?iso distribution roughly matches our 
simulated distribution. 

The L^pk.obs distribution for all simulated bursts (solid 
line) and detected bursts (dashed line) is shown in Figure 

t Similar to the case of the detected i?iso distribution, 
e low end of the detected i?pk,obs distribution is largely 
influenced by the limited sensitivity of the detector below 
^ 25 keV, whereas the high end of the detected Spk.obs 
distribution reflects the true cut-off in the simulated pop- 
ulation. 

In Figure [Sj I plot E'pk.obs versus the observed energy 
fluence S'boi, integrated over an energy range of 10 to 
10000 keV, for our entire simulated data set (dots), along 
with the bursts that would have triggered the instru- 
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ment (circles). The color of each data point represents 
the simulated burst's trigger significance in the 50-300 
KeV energy range as seen by the detector, which can be 
seen to correlate with the observed energy fluence. The 
boundary between the detected and undetected popula- 
tions tracks the energy dependance of the detector's ef- 
fective area, with the weakest detections occurring near 
the center of the detector's energy window, where the 
instruments sensitivity is maximized. By following the 
color gradient, it can be seen that increasing the detec- 
tion threshold acts to narrow the observed population. 
The results presented in FigurelSlare in qualitative agree- 
ment with similar modeling of the BATSE t hresho ld as 
a function o f energy found by Nava et al.| (|2008|) 



As 

discussed by |Nava et al. (2008), 1 hnd that the lack of 
detected GRBs m the high 6'boi, low i?pk,obs regime of 
Figure |8] is not due to any limitations of the detector, 
but rather starts to reflect the intrinsic cutoff of the GRB 
population as we run out of bright, nearby events. 

To understand the nature of this cutoff in the energy 
fluence distribution, we turn to Figure|9) where I plot the 
source frame Ep]^ versus the isotropic equivalent photon 
luminosity Lisa for our entire simulated data set (dots), 
along with the bursts would have triggered the instru- 
ment (circles). The color of each data point now rep- 
resents the simulated burst's redshift. A broad scatter 
plot seen when considering the entire data set reflects the 
fact that I assume no intrinsic correlation between Liso, 
Spk, or redshift in the process of generating the simulated 
population. The sharp drop in the number of events at 
the right end of the plot reflects the relatively steep lu- 
minosity function assumed in our population synthesis 
code, resulting in the cutoff in the S'boi distribution seen 
in Figure [8] 

Several patterns become apparent when I consider only 
the detectable events, which form a triangular region in 
the i?pk versus Liso plot. Although no correlation is 
present between Liso, ^'pk, and redshift for the entire 
sample as a whole, a very distinct pattern emerges when 
considering only the detected events. These patterns 
are due to the energy dependance of the instrument's 
detection threshold. The decreasing effective area near 
the high and low bounds of the instrument energy win- 
dow becomes increasingly important for low luminosity 
bursts, decreasing the detection efficiency for low lumi- 
nosity bursts with very low and very high i?pk,obs values. 
On the other hand, extremely bright bursts are relatively 
easier to detect, even if -Epk.obs is outside of the detec- 
tor's energy window, hence the wider range of detected 
bursts on the right side of Figure [9j Finally, although 
intermediate luminosity bursts are seen at all redshifts, 
we only see the most luminous events at high redshift, 
again due to the Malqucst bias. More importantly, lu- 
minous bursts with low _Epk are not seen at very high 
redshift since their -Epk.obs values would be redshifted 
well outside of the detector's energy window, resulting 
in a low observed flux and making them less likely to be 
detected. These two factors, a Malquest bias and the en- 
ergy dependent detection threshold, produces a redshift 
gradient th at is correlated with both Ep\. and Liso- 

In Figure [Toj I plot Eiso versus -Epk for our entire sim- 
ulated data set (dots), along with the bursts that would 
have triggered the instrument (circles) . The color of each 
data point again represents the burst's redshift. Again, 
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Fig. 8. — The observer frame i?pk,obs versus the observed energy 
fluence Sboi for simulated (dots) and detected (circles) bursts. The 
color gradient represents trigger significance of each burst in the 
50-300 KeV energy range. Increasing the detection threshold acts 
to narrow the observed population. 
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Fig. 9. — Source frame Bpk versus the isotropic equivalent pho- 
ton luminosity Liso for all simulated (dots) and detected (circles) 
bursts. The color of each data point represents the simulated 
burst's redshift. Although no correlation is present between Liaoi 
-Epk, and redshift for the entire sample as a whole, a very distinct 
pattern emerges when considering only the detected events. 

a broad scatter plot between the simulated Eiso and i?pk 
values for the entire data set (detected and undetected) 
reflects the fact that I assumed no intrinsic correlation 
between a GRB's Spk, photon luminosity iiso, and du- 
ration in the process of generating the simulated pop- 
ulation. Despite this, a broad positive correlation does 
appears between Eiso and -Epk among the GRBs that 
would have been detected by the instrument. 

The slope of the lower-right edge of this correlation 
is roughly linear and can be understood as occurring by 
definition of the Band spectral model used in our simula- 
tions. The Band function is comprised of two power-law 
components smoothly joined by an exponential. iJpk is 
defined as the location of the break between these com- 
ponents and reflects the peak of the vF^, spectra. So, 
for GRBs with identical Liso, the events with higher _Epk 
values will by definition exhibit larger energy luminosi- 
ties and hence larger Eiso, modulo the GRB duration. 
This can be seen more easily in Figure where I plot 
-Epk versus i?iso with a color gradient representing Liso- 
Areas of constant Liso can be seen to track a linear corre- 
lation between _Epk and Eiso- So the lower-right edge that 
appears in Figure [lO] refiects the sharp cutoff in our pre- 
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Fig. 10. — Source frame iJpjj versus Eisa for all simulated (dots) 
and detected (circles) bursts. The color of each data point repre- 
sents the simulated burst's redshift. A broad positive correlation 
appears between E[g^ and Spj; when considering only the detected 
population. 
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Fig. 11. — Source frame Epi^ versus Eiao with a color gradient 
representing Liso- Areas of constant L^^^ can be seen to track a 
linear correlation between Bpi^ and Ei^o due to the conversion from 
photon luminosity to total energy, modulo the burst duration. 

scribed luminosity function above Lcut ~ 53 transformed 
into an energy luminosity, which is inherently correlated 
with the burst's spectra. Likewise, the upper-left edge of 
the observed correlation reflects the shape of the detec- 
tor threshold observed between Liso and ii^pk in Figure 
[9| but now transformed from photon luminosity to total 
energy, again modulo the GRB duration. 

As with the photon luminosity, there is a very distinct 
correlation between redshift, Eiso, and £'pk, with increas- 
ingly energetic and hard events only being detected at 
high redshift. Burst's with low Eiso are only detectable 
at low redshift, with only a slight dependance on Epy^ as 
described above. Bursts with both high Eiso and high 
£^pk are rare and therefore the GRB luminosity function 
dictates that these events are most likely to occur at 
higher redshift, due to the increase in the observed vol- 
ume. Bursts with high Eiso and low Ep\^ are also rare, but 
become increasingly difficult to observe at high redshift 
because their iJpk.obs will be redshift below the detector's 
energy window. In fact, low iJpk of all energies suffer the 
same selection bias, relegating low E'pk events to low and 
intermediate redshifts. The net effect of these selection 
biases is a modification of the slope correlating Eiso and 
£'pk from the linear relationship expected from the na- 



ture of the Band function to a shallower slope due to the 
depletion of burst's with intermediate to high Eiso and 
low E'pk values. 

We can see this effect more easily in the left panel 
of Figure [12] where I again plot E'iso versus £^pk for our 
entire sirnulated data set (dots), along with all bursts 
that would have triggered the instrument (circles). The 
color coding now represents the trigger significance in the 
50-300 KeV for the detected population. Although there 
are weakly detected GRBs throughout the detected range 
of i^iso versus Ep\., the bursts with the highest trigger 
significance start to form a very distinct pattern. The 
edge of the luminosity function becomes less prominent 
through the depletion of low i?pk values with decreasing 
ii'iso, with the majority of bursts cluster at intermediate 
ii^iso and £^pk values, producing a slope which is shallower 
than the linear relation between Eiso, and iJpk for a given 

.^iso ■ 

I quantify the net effect these selection biases have 
on the observed GRB population by converting the left 
panel of Figure [12] into the 2-dimensional probability dis- 
tribution. By mapping the intensity of this probability 
distribution, I can determine the GRB population most 
likely to be detected by the instrument. The resulting 
probability distribution can be seen in the right panel of 
Figure [12] 

The intensity gradient and associated contours repre- 
sent the most likely set of Eiso and i?pk values to be 
detected by the instrument once all of the various selec- 
tion biases discussed above have been taken into account. 
For comparison, I have plotted the slope, normaliz ation 
and scatt e r of th e £^pk — i?iso correlation reported by Am 



ati et al. (2008). The slope of the reported correlation 



is roughly a ~ 0.54, matching the slope of the intensity 
plot produced from our simulated data set. The slope 
of the probability distribution is governed by the edge 
of the luminosity function at high ii^iso and £'pk, but de- 
viates at low £'pk due to the depletion of these events 
because of cosmological redshift. 

5. DISCUSSION 

The GRB model presented here is a very simplified 
representation of the complex nature of GRBs. All of 
the bursts produced by our simulation are single pulsed 
events with Band function representations of their pho- 
ton spectra. Furthermore, the code explicitly assumes 
that their spectral evolution is dictated by simple rela- 
tivistic kinematics with a pronounced hard to soft evo- 
lution. In reality, GRBs are typically described by far 
more complex time histories, exhibit a wider range of 
spectral shapes, and at times do not follow the hard to 
soft evolution predicted by relativistic kinematics. Still, 
these simulations serve to demonstrate how selection ef- 
fects caused by a combination of instrumental sensitivity 
and the cosmological nature of GRBs can act to produce 
a strong correlation between observed properties. 

As with all flux limited observations, our knowledge 
of astrophysical sources will be heavily influence by the 
properties of the brightest events in the sample, since 
they are the easiest to detect, even if they are not the 
most common. T he r esults presented in the left and right 
panels of Figure ]T2] show that this is no different for 
GRBs, with the brightest bursts in the GRB luminos- 
ity function being the most commonly detected events. I 
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Fig. 12. — Left Panel: Ei^o versus i?p]j for all simulated (dots) and detected (circles) bursts. The color of each point represents the burst's 
trigger significance in the 50-300 KeV for the detected population. The bursts with the highest trigger significance start to form a very 
distinct pattern. ijight Panel: The conversion of the left panel into a 2-dimensional probability distribution representing the most likely set 
of i?iao and Bpi; values to be detected by the instrument once all of the various selection biases have been taken into account. The dashed 
and dotted line represent the slope, normalization, and 1 sigma scatter of the Ep]^ — E^so correlation reported by |Amati et al.| ( |2008[ l. 



find that the edge of the GRB luminosity function plays 
an important role in producing the boundaries of the 
observed correlation between i?pk and E'iso- The fact 
that no instrument has to date detected a GRB with an 
-E-iso > lO'^^ erg is due to a cutoff in the source population 
and cannot be explained by any selection effects. There- 
fore, the upper right region of the observed Epy^ — Eiso 
diagram reflects the true cutoff in the GRB population. 
Despite this, our simulations show that this cutoff, when 
expressed in energy space, cannot alone explain the slope 
of the -Epk — i^iso relation. Instead, the non-detection of 
low Epk, high Eiso GRBs at high redshifts acts to flatten 
the slope of the relationship. By examining the density 
of detected bursts in the i?pk — -Eiso parameter space, the 
simulations reveal that the area in which GRBs are most 
likely to be detected roughly ni atches the a ~ 1/2 slope 



(2008) 



reported by Amati et al 

Still there are several facets of the E'pk — i^iso that 
are not addressed by these simulations. Foremost, I 
do not address the nature of X-ray flashes (XRFs) and 
the adherence of at least three XRFs with i^pk obs < 20 
KeV (XRF 020903, XRF 050416, and XRF 060218) to 
the correlation. Events with such low £'pk.obs values 
were not simulated in our analysis because of BATSE's 
lack of energy coverage at those energies. Furthermore, 
the instruments that detected these events, BeppoSAX 
for XRF 020903 and Swift-BAT for XRF 050416, and 
XRF 060218, have unique response functions and trig- 
gering criteria which prevent me from quantifying the 
observed parameter space in which these events are most 
likely to occur. 

Independent of detector consideration though, if XRFs 
are drawn from the same photon luminosity function as 
their higher energy G RB counterparts, as sug gested by 
Lamb et al. (20031 and Sakamoto et al. (2008), then low 
i/pk events would by dehnition have low\fc/'iso for a given 
photon luminosity, modulo the burst duration. An XRF 
with i?pk ^ a few keV will not posses the same energy 
luminosity as a burst with E'pk ~ 500 KeV of the same 
photon luminosity. Therefore, the cutoff in the photon 
luminosity function that is responsible for the production 
of the edge in the i?pk — i^iso diagram would, by defini- 
tion, preclude the existence of high E^so, fow iJpk events 
unless XRFs are drawn from a separate population with 



a very different photon luminosity function compared to 
cosmological GRBs. 

If XRFs, XRRs, and GRBs do represent a continuum 
of events drawn from the same population, then Figure 
[8] would indicate that soft and weak XRFs should have 
an equal number of hard and weak counterparts with 
Eu k > 100 keV. Such eve nts are not in fact reported 
by Sakamoto et al. (20081, where the authors compile 



a comprehensive list of XRFs and XRRs detected by 
Swift-BAT down to Se ~ 10""^ erg cm'^ in the 15-150 
keV energy range. Hard events at these fluence levels 



are reported, though, by Butler et al. (2007), who utilize 
a novel Bayesian spectral fitting technique to estimate 
£'pk.obs near or beyond the upper edge of the Swift-BAT 
energy window. Therefore, the results presented in Fig- 
ure[8]which point to the existence of hard and soft events 
at low fluence levels are largely consistent with the re- 
sults found by Butler et al. ( 2007 ) . If GRBs have a broad 
photon luminosity function and iJpk distribution and oc- 
cur at a wide range of redshifts, then it is expected that 
the vast majority of bursts seen in the observer frame 
should occupy the low ^boi regime while covering a wide 
range of Epy^ values. 

Only a small fraction of the GRBs that form the basis 
for the i?pk — -Eiso correlation were detected by BATSE, 
but the instrument's broad energy range and its rela- 
tively straight forward trigger criteria make it extremely 
easy to model for this analysis. Furthermore, the pur- 
poses of the analysis presented here is not necessarily 
to match the exact form of the -Epk — -Eiso correlation 
observed by other instruments, but rather to demon- 
strate how the combination of detector thresholds and 
the population distribution can work to produce artifl- 
cially strong correlations that neither exist in the source 
frame of the population or are tied to the underlying 
physics of the sources in question. 

6. CONCLUSIONS 

The simulations present here paint a complex picture 
regarding the nature of the observed correlation between 
-Eiso — Epk. By simulating a population of GRBs using an 
assumed luminosity function and co-moving rate density, 
I find that 

• The Epk, Eiso, and Ejso distributions among the 
detected population can appear strongly correlated to 
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the observer despite the existence of only a weak and 
broad correlation in the original simulated population. 

• The energy depcndance of any detector's flux limited 
detection threshold acts to produce a correlation be- 
tween the source frame i?pk and Ei^o for low luminosity 
GRBs, producing the left boundary of the observed 
Epk — Eiso correlation. 

• Very luminous GRBs arc found at higher redshifts 
than their low luminosity counterparts due to the 
standard Malquest bias, causing bursts in the low £^pk, 
high Eiso regime to go undetected because their iJpk 
values would be redshifted to energies at which most 
gamma-ray detectors become less sensitive, producing 
the right boundary of the observed correlation and 
flatten the slope of the relationship. 

• The origin of the observed correlation is a complex 
combination of the instrument's detection threshold, the 
intrinsic cutoff in the GRB luminosity function, and the 
broad range of redshifts over which GRBs are detected. 



Although the GRB model presented here is a very sim- 
plified representation of the complex nature of GRBs, 
these simulations serve to demonstrate how selection ef- 
fects caused by a combination of instrumental sensitivity 
and the cosmological nature of an astrophysical popula- 
tion can act to produce an artificially strong correlation 
between observed properties that does not exist in the 
source frame of the population. 



D.K. acknowledges financial support through the 

Fermi Guest Investigator program. We thank Nathaniel 
Butler. Robert Prcecc. and Josh Bloom for insightful dis- 
cussions. This work was supported in part by the NASA 
Fermi Guest Investigator program under NASA grant 
number NNX10AD13G and the U.S. Department of En- 
ergy contract to the SLAG National Accelerator Labora- 
tory under no. DE-AC3-76SF00515. 



REFERENCES 



Amati, L., Prontcra, F., Tavani, M., in't Zand, J. J. M., 

Antonelli, A., Costa, E., Feroci, M., Guidorzi, C, Heise, J., 
Masctti, N., Montanari, E., Nicastro, L., Palazzi, E., Pian, E., 
Piro, L., & Soffitta, P. 2002, A&A, 390, 81 
Amati, L., Guidorzi, C, Prontcra, P., Delia Vallc, M., Pinelli, P., 

Landi, R., & Montanari, E. 2008, MNRAS, 391, 577 
Band, D., Mattcson, J., Pord, L., Schacfcr, B., Palmer, D., 

Teegarden, B., Cline, T., Briggs, M., Paciesas, W., Pendleton, 
G., Pishman, G., Kouveliotou, C, Meegan, C, Wilson, R., & 
Lestrade, P. 1993, ApJ, 413, 281 
Band, D. L., & Preece, R. D. 2005, ApJ, 627, 319 
Butler, N. R., Bloom, J. S., & Poznanski, D. 2010, ApJ, 711, 495 
Butler, N. R., Kocevski, D., & Bloom, J. S. 2009, ApJ, 694, 76 
Butler, N. R., Kocevski, D., Bloom, J. S., & Curtis, J. L. 2007, 
ApJ, 671, 656 

Pishman, G. J., Meegan, C. A., Wilson, R. B., Brock, M. N., 
Horack, J. M., Kouveliotou, C, Howard, S., Paciesas, W. S., 
Briggs, M. S., Pendleton, G. N., Koshut, T. M., Mallozzi, R. S., 
Stollberg, M., & Lestrade, J. P. 1994, ApJS, 92, 229 



Gehrels, N., Cliincarini, C, Gionmii, P., Mason, K. O., Nousek, 
J. A., Wells, A. A., White, N. E., Barthelmy, S. D., Burrows, 
D. N., Cominsky, L. R., Hurley, K. C, Marshall, P. E., 
Meszaros, P., Roming, P. W. A., Angelini, L., Barbier, L. M., 
Belloni, T., Campana, S., Caraveo, P. A., Chester, M. M., 
Citterio, C, Cline, T. L., Cropper, M. S., Cummings, J. R., 
Dean, A. J., Peigelson, E. D., Penimore, E. E., Prail, D. A., 
Fruchtcr, A. S., Garmirc, G. P., Gcndreau, K., Ghisellini, G., 
Grcincr, J., Hill, J. E., Hunsbcrgcr, S. D., Krimm, H. A., 
Kulkarni, S. R., Kumar, P., Lcbrun, P., Lloyd-Ronning, N. M., 
Markwardt, C. B., Mattson, B. J., Mushotzky, R. P., Norris, 
J. P., Osborne, J., Paczynski, B., Palmer, D. M., Park, H.-S., 
Parsons, A. M., Paul, J., Recs, M. J., Reynolds, C. S., Rhoads, 
J. E., Sasseen, T. P., Schacfcr, B. E., Short, A. T., Smale, 
A. P., Smith, I. A., Stella, L., TagUafcrri, G., Takahashi, T., 
Tashiro, M., Townsley, L. K., TucUer, J., Turner, M. J. L., 
Victri, M., Voges, W., Ward, M. J., Willingale, R., Zerbi, 
P. M., & Zhang, W. W. 2004, ApJ, 611, 1005 

Ghirlanda, G., Nava, L., Ghisellini, G., Firmani, C, & Cabrera, 
J. I. 2008, MNRAS, 387, 319 

Kocevski, D., Ryde, F., & Liang, E. 2003, ApJ, 596, 389 

Lamb, D. Q., Donaghy, T. Q., & Graziani, C. 2003, ArXiv 
Astrophysics e-prints 

Meegan, C, Lichti, G., Bhat, P. N., Bissaldi, E., Briggs, M. S., 
Connaughton, V., Diehl, R., Pishman, G., Greiner, J., Hoover, 
A. S., van der Horst, A. J., von Kienlin, A., Kippen, R. M., 
Kouveliotou, C, McBreen, S., Paciesas, W. S., Preece, R., 
Steinle, H., Wallace, M. S., Wilson, R. B., & Wilson-Hodge, C. 
2009, ApJ, 702, 791 

Nakar, E., & Piran, T. 2005, MNRAS, 360, L73 

Nava, L., Ghirlanda, G., Ghisellini, G., & Firmani, C. 2008, 
MNRAS, 391, 639 

Preece, R. D., Briggs, M. S., Mallozzi, R. S., Pendleton, G. N., 
Paciesas, W. S., & Band, D. L. 2000, ApJS, 126, 19 

Sakamoto, T., HuUinger, D., Sato, G., Yamazaki, R., Barbier, L., 
Barthelmy, S. D., Cummings, J. R., Penimore, E. E., Gehrels, 
N., Krimm, H. A., Lamb, D. Q., Markwardt, C. B., Osborne, 
J. P., Palmer, D. M., Parsons, A. M., Stamatikos, M., & 
Tueller, J. 2008, ApJ, 679, 570 



